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1. INTRODUCTION 

The three-phase power distribution system is one of the most popular in the world, and it allows to 
transfer of energy from production to consumption much more efficiently than single-phase power 
transmission systems. One of the problems of three-phase systems is that the power consumption is not 
evenly distributed between the phases, as a result of which the system becomes unbalanced, causing an 
increase in energy losses and damage to electrical equipment. The imbalance in the three-phase systems is 
one of the most critical issues of the three-phase system, it has been extensively studied for many years, and 
specific achievements have been achieved in solving this problem [1]-[5]. Thanks to computerization, it has 
become possible to apply highly sophisticated and automated methods to improve the quality of the systems. 
The negative sequence currents create a magnetic field that rotates (against the direction of rotation of the 
synchronous machine), which retards the rotor. These fields also induce currents of doubled frequency in the 
rotor of machines with a cylindrical rotor at the poles of the salient-pole machines. These currents induced 
overheat the rotor. Heating, in turn, leads to a loss of mechanical integrity or failure of the electrical 
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machine's insulation. In asynchronous motors, even a 5% imbalance can result in a 25% reduction in engine 
power. The imbalance presents in the supply voltage by 3% can increase the heating of the rotor by about 
20%, so three-phase machines must be protected against unbalanced currents [6]. 

For some applications, it is possible to reduce the imbalance by changing the operating parameters. 
To reduce the effect of negative sequence currents, a cause of the negative sequence voltage drops, the 
supply voltage requires a low impedance of the internal system. However, applying this method in rural 
networks is complicated because it is rather challenging to select a single-phase, automatically switched the 
load of the appropriate power. Therefore, it is not possible to significantly reduce the load asymmetry. In 
addition, when using balanced devices of this type, the reliability of the power supply of single-phase 
switchable load decreases; The measuring and commutation complex of these devices seems rather 
complicated. Another mitigation method is the use of special transformers such as a three-phase or three- 
phase transformer or Scott and Steinmetz transformers [7]-[10]. Transformer balanced three-phase can 
reduce energy losses by reducing the amplitudes of harmonics, reducing the resistance. This increases the 
working life of energy sources in networks with phase misalignments. The devices are designed to increase 
the reliability of autonomous generators and consumers when the loads are asymmetrical and for balancing 
phase voltages by decreasing the resistance of the zero sequence, so that the same phase bias in the voltage is 
eliminated, which happens mainly due to the asymmetry of the load currents phases. Reactive power 
compensation can be performed using a synchronous generator to generate or collect reactive power, but it 
has that disadvantage: the high cost of installation, maintenance, and complex operation. Synchronous 
generators are often used to concentrate compensation with a large capacitance, and a compensation 
capacitor overcomes the disadvantages of a synchronous generator, but the drawback is only for use with the 
system's bandwidth small, sensitive to voltage changes. 

Combined balancing devices always give the best result, however, they are becoming more complex 
in the calculation, and the number of compensating components necessary for total compensation requires 
more. The cost of the compensation system also increases. Considering the growing demand for power 
systems for measures to improve energy quality and reduce the cost of investments and maintenance costs of 
power quality improvement systems, we need proposals that ensure low-cost symmetrization of three-phase 
systems [10]-[15]. The complexity of consumption patterns causes difficulty and even impossibility to create 
a general method for optimizing the unbalanced system. The design of universal symmetrizing systems 
requires a combination of several approaches and several devices. This leads to an increase in the cost of 
creation and maintenance. For example, in many cases, for industrial load, the repeatability of load 
characteristics (including their asymmetry) is observed at different time scales from daily to annual. This 
circumstance makes it possible to simplify and reduce the cost of symmetrization systems, which is essential 
since investment costs are the user's primary concern within the acceptable quality of electric systems. With 
these analyses, we describe the methodology for optimizing three-phase systems operating in asymmetric 
modes. In meeting the requirements for the quality of the three-phase power system and the investment costs, 
this approach aims to reduce the number of reactive power compensation devices [15]-[20]. 


2. RESEARCH METHOD 
2.1. New proposed fractional-polynomial functions 

Assume that we have a three-phase circuit consisting of (n+1) nodes and m (n+1<m) branch and 
described by the matrix of link nodes [10-20]: 


Q11 2 Qim 
a21 Q22 A2m 


Anı An2 Anm 
Vector of the conductance of the branches is 
Y = diag, Yo, ... Ym) 
Vectors of current and electromotive force sources is formulated as 
J = UJ Smd 
E = (E, E, ... Em)t 


The node voltage equations are written in matrix form as: [18] 
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AYA'U, = —A(J + YE) (1) 
The node voltages can be found by the Cramer formula [19]-[21]: 


= det Y; 
tT det Y 


The matrix Y will be: 


1 Yui Yi Vij Yin 
1 Yia Ya Yay Yin 
J Yj Yi Yy Vin 
n Yn Yni Ynj Yan 
The matrix determinants of Y and Y; are defined as follows: (here and were with sign (*), are results 


of the circuit analysis we have done on the maple software, and due to the limitations of paper we will not 
present in more detail). 


ay a2 a3 
det Y = ao + A,X, + a2X3 + a3X1X2 = Ay (1 +— x, +— xX + x2) 
ao ao ao 
det Y; = boi + bi iX4 + bz iX2 + b3 iX1X2 
In general, voltages will be: 


a det Y; boi + Dy iX1 + by iX2 + b3 iX4X2 


= = 
det Y ( ay a2 a3 ) 
ao I+” ta” tga 


Let a, = ap/a0; p = 1 +3 and c4, = bgi/do3; q = 0 +3, we have: 


det Y; CoitCy iX1+C2 {X2 +C3 įiX1X2 
U; = i — vitei it it (2) 


det Y 14+04X14+02X%24+03X4X2 
If the i-th branch connects two nodes j-th and k-th, the current flow in that from k to j, in general, as: 


p= do i+d1,iX1+d2,iX2+d3 iX1X2 (3) 
d 1+a1X1 +A2X2+QA3X1X2 


Coefficients a, ...3 in this case has the same value as the coefficients of the voltages, i = 1 + m — 
index of branches [20]-[24]. 


2.2. Testing 
In the second experiment (proposed functions), we performed as described in section 2. The 
example circuit is shown in Figure 1. The error between two experiments was calculated by the formula: 


|o @1%2)|-[9, @1x2)| 


e (x1,x,) = 100% 
9, CeLx2) 


here gs Cx x2) = jee E or Oe? Cex) bik = 0,4;j + k;m = 1,2 — currents (voltages) of 
the first experiment and the proposed method [25]-[33]. 

In Figure 2 and Figure 3, present the errors €,;(%,, x2) of current and €y(x,, x2) of the voltage of the 
proposed FPF without considering random errors (€, = 0), such as errors of measurement or constant change 
of load etc. However, this case does not exist in practice, so we had taken into account random error in the 
range of £, = +0.5% and +1%. The errors were presented as in Figures 4 to 7. The root means square 
(RMS) values of the errors were calculated by the formula, and the results will be contained in the Table 1. 


11 O_o 
_ iise ya 22 
Erns = A|>, pene (x! Xo 


here p,q - the lengths of vector x1, x3. 
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FACTS! 
"gh Z4=1j 

oS 5 

Z1 = 0.01 Ohm (4) 
E= 1.e(4.2.pi/3) 27 = 0.01) 

Q) A =] 

Z Z2=0.01 Ohm 75 = 142) 
Ec=1.e().2.pi/3) G) 
O T = — 

Z3 = 0.01 Ohm z014 


Z10 = 242) 


FACTS2 


Figure 1. Example circuit 


Error of 11 (%) 
Error of U1 [%] 


x2 [Ohm] 5 2 ct form x2 [Ohm] K x [Ohm] 
Figure 2. Dependence of error in determining the Figure 3. Dependence of error in determining the 
’ 
value of the error £; on reactance’s values x4, X2, value of the error €y on reactance’s values x4, x2, 
= 09 
when £, = 0% whens % 


Error of 1 (%) 
Error of U1 [%] 


~ 2 as — x 2 
x2 [Ohm] oo x1 [Ohm] x2 [Ohm] ° oO x1 [Ohm] 
Figure 4. Dependence of error in determining the Figure 5. Dependence of error in determining the 
value of the error €; on reactance’s values x1, X3, value of the error £y on reactance’s values x4, X3, 
when €, = 0.5%. when £, = 0.5 %. 
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Error of 1 (%) 


aT 7 x2 (Ohm) 
x2 [Ohm] 


a 2 
o o x1 [Ohm] 


Figure 6. Dependence of error in determining the Figure 7. Dependence of error in determining the 
value of the error £; on reactance’s values x4, X3, value of the error £y on reactance’s values x4, X2, 
when €, = 1% when £, = 1% 


Table 1. The RMS values of the errors depend on €, 
Er, % 
0 0.5 1 
E rus, % 6.87-10% 0.34 1.14 
Ermsy% 9.02-10 ° 0.73 7.98 


We had repeatedly executed it (n =1000 times) to assess errors and accurately received average errors. 


E (x2) = EZ E(x x2) [%] 


n 


The RMS values of the average errors are contained in Table 2. 


where p, q - the lengths of vector x4, x2. 


Table 2. The dependencies of the root mean square errors on €,. 
En% 
0 0.5 1 
E pus, % 6.87- 10° 1.31 2.59 
Ermsu%  9.02-10° 158 3.01 


As we can see, the most significant error of the proposed method is under 2%, when the most 


significant random error of £, = 0.5% and about 3% when £, = 1%. The accuracy of fractional-polynomial 
functions depends on the error of measurement. Nowadays, the development of measurement devices with 


high quality makes a tiny error €, < 0.5%, however, in this article, we not only take into account the errors 
of measurement but also including the constant change of load as shown in Figure 8. 


"dZ < DZ 


az Z(ti) 


Figure 8. The constant change of loads 
Permissible change limits of load: 
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g, = ZED- _ DZ 100% 


k Z(to) Zo 
dZ = |Z — Zol 


where Zo = Z(to)- impedance at the beginning of the calculation; Z(t,)-impedance at the end of the 
calculation. At moment t4, the calculation process is completed, and the values (x1, X2) received will be 
applied to the compensation system if dZ < DZ. 


3. NUMERICAL RESULTS AND DISCUSSION 
Sequence components are formulated as [20]-[37]: 


[Ok 
1 
i Cx, x2) TO C, x2). (4) 


(j,k 
) (x1, X2) 1 1 177 i9 (x1, x2) 
(i) a @ 1 Gk) 
Ip (x1, x2) Ig (x1, x2) 


where j,k = 1, N; j + k; We got: 


j,k j,k j,k j,k 
ae ene x +c xatd YO X41 X2 


:(j,k 
i9 (x1, x2) = 


(5) 


14+04%X%44+02%2+03X14X2 


Negative and zero sequence components appear in asymmetrical systems, and they have greater 
values when the system is more asymmetric. 


» Uk) » Gk) . 
|L : (x1, X2)| > max f= 1/th 7 (x1, x2)| > min 
ILI Gen xa) > minor 4 fy = hI” Œ) > min © 
+ Gk) ‘ + Uk) ; 
llo (x1, X2)| > min fz = lo (x1, X2)| > min 


Let F(X) = {fa fo, fa} 


The multi-objective optimization problem can be brought back to mathematical form and described 
as follow [2], [12] 


min F (x1, x2) = (fi x1 x2), f2 1, X2), fa (X1, X2)} d 
s.t. (%4,X2) E X, 


where the integer k > 2 is the number of objectives, and the set X is the feasible set of decision vectors. The 
valued vector of the objective function is defined as: 


F:X > R* 


Therefore, X_1 is anon-dominated solution and X_2is an inferior solution. The solution X_1 will be 
called dominated solution if: 


Wi = 1..3: f; (X1) < f; (X2) 


Typically, in the multi-objective optimization problem, we often encounter non-denominated 
solutions rather than denominated solutions. As we can see in Figure 9 in 2-D coordinate system. The most 
important relationship between all objectives is the Pareto front, and we can choose the appropriate solution 
as the optimal solution of the problem when we gave a weight coefficient to each objective [22]-[27]. For 
multi-objective optimization problems solved by genetic algorithms. The set of solutions after every single 
loop of the iteration is called generation. An initial set of solutions is randomly selected and designated as the 
first generation. In every iteration loop, the genetic operators, selection, crossover, and mutation are applied 
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to the previous solutions (previous generation) to generate the new solutions. The final set of solutions, also 
called Pareto font, was obtained after a certain number of iterations [15]. 


10 


\—/T1 110/10 
(T) kv 


—1 JT 7 L 


(A) (A) (A) (A) (A) 
\ST1 100.4 \S/T2 100.4 T3 10/04 \S/T4 100.4 \S/TS 10/0.4 
(K) kV ) kV KA kV A) kV (A) kV 
— — — — — 
1 4-Lines 4-Lines 4-Lines 4-Lines 4-Lines 
> — — == — — 
o = = = == PF =m A FACTS 1,2 FACTS 3,4 FACTS 5,6 FACTS 7,8 FACTS 9,10 
m 10 Load 1 Load 2 Load 3 Load 4 Load 5 
Figure 9. Pareto front in 2-D coordinate system Figure 10. Electrical system of the glass factory 


4. 


CONCLUSION 
In this paper, a new method to overcome these problems is proposed and investigagted. The 


proposed method is based on the fundamental electrical quantities (voltages and currents) on the controllable 
values of the static compensation devices and overcoming of the asymmetric operation regime in the three- 
phase systems. Furthermore, the efficiency and method speed are improved by intelligent search algorithms 
(NSGA, GA). The research results show that the new mothed is powerful and efficiency. 
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